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We study reconstruction of the lensing potential power spectrum from CMB temperature data, 
with an eye to the Planck experiment. We work with the optimal quadratic estimator of Okamoto 
and Hu, which we characterize thoroughly in application to reconstruction of the lensing power 
spectrum. We find that at multipoles L < 250 our current understanding of this estimator is biased 
at the 15% level by beyond-gradient terms in the Taylor expansion of lensing effects. We present 
the full lensed trispectrum to fourth order in the lensing potential to explain this effect. We show 
that the low-L bias, as well as a previously known bias at high-L, is relevant to the determination 
of cosmology and must be corrected for in order to avoid significant parameter errors. We also 
investigate the covariance of the reconstructed power, finding broad correlations of ~ 0.1%. Finally, 
we discuss several small improvements which may be made to the optimal estimator to mitigate 
these problems. 



I. INTRODUCTION 

The cosmic microwave background (CMB) provides us with a picture of the Universe at the time when neutral 
atoms first formed, some 400, 000 years after the big bang. The O(10~^) anisotropies in its temperature characterize 
the density perturbations at that time, and from them we have been able to learn a good deal about the properties of 
our Universe (see Ref. [T] for the most up-to-date constraints). The CMB is a powerful tool because it is largely free 
of astrophysical complications. Its features may be derived from initial conditions using linear physics, and to a good 
approximation it may be treated as a Gaussian random field statistically characterized by a single power spectrum. 

At a high level of sensitivity and resolution, however, one expects to find small deviations from Gaussianity. These 
are expected at some level both from the high-energy physics which generated the primordial perturbation thus seeding 
the CMB anisotropies, and from secondary astrophysical effects after the CMB was released at recombination (see 
Ref. [2 for a review). One of the most important astrophysical non-Gaussianities is that introduced by gravitational 
lensing of the CMB (see Ref. [3] for a review.) We view the distant Universe through the intervening matter of our 
Hubble volume, the gravitational potential of which distorts the paths of the CMB photons. In essence, we view the 
CMB through bubbled glass. This blurs the features of the underlying anisotropies and introduces non-Ganssianity 
into the observed sky [5j. In some sense, lensing of the CMB is a nuisance as it provides a potential contaminant for 
measurements of intrinsic non-Ganssianity in the early Universe O [7] , with which we hope to be able to discriminate 
between fundamental theories. From another perspective, however, CMB lensing is a tool which provides us with a 
new perspective on the contents of our Universe ^ (9j . By studying the warping effect which large scale structure has 
on the paths of CMB photons, we may make integrated measurements of the matter power spectrum at intermediate 
redshifts. This provides us with additional powers to refine parameter measurements, and to break degeneracies in 
measurements of the primary CMB [TUHH]. With Planck lensing data, for example, we expect to be able to constrain 
the sum of neutrino masses to sub-eV levels [13] . 

The effect of lensing is predominantly a simple remapping of the temperature distribution of the underlying, unlensed 
CMB by the lensing deflection. To leading order, the deflection is the gradient of a lensing potential <j> which is an 
integrated measure of the intervening gravitational potential along the unlensed line of sight (e.g. [3j)- We may write 
the remapping on the sphere as 

e(fi) -e[n + V0(n)], (1) 

where G>(n) is the unlensed CMB temperature in the direction n and Q{n) is the observed, lensed temperature.^ 
On the scales which we shall be considering, (/)(n) is well approximated as a Gaussian random field, and can be 
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^ Equation (nl is rather symbolic; the remapping is by a distance |V(/)| along the geodesic tangent to V</> through n |14| . 



2 



characterized solely by its power spectrum Cf/' in harmonic space. This power spectrum is calculable for a given 
cosmology (see [3] for details) and can be compared to measurements in order to constrain parameters. 
We can gain insight into the effects of lensing by considering the Taylor expansion of Eq. ([T]): 

e(n) = e(n) + V,0(n)V^e(n) + --- , (2) 

from which we can see that the effect of lensing is to introduce dependence on the higher derivatives of the unlcnsed 
CMB temperature into the CMB itself. These additional terms can couple to give non-Gaussian effects on the sky. To 
the extent that the lensing deflections are small, the first order term in Eq. ^ is sufficient to describe the lensing (but 
see |15j for a discussion of circumstances where this approximation fails). Therefore, one expects that an estimator 
of the form 

V0(n) oc e(n)Ve(n) (3) 

will give a reasonable reconstruction of the deflection field. Taking the expectation value over CMB realizations we 
can see that if properly normalized this estimator recovers the deflection field in the mean (at first order in the lensing 
expansion). The reconstruction is clearly quite noisy, however, as for any particular realization there will be some 
spurious 0V0 correlation. This term is zero in the mean, but contributes significant variance to the estimator. A 
further issue is what field to use as a proxy for the unlensed CMB. 

We can achieve a somewhat better overall reconstruction than that in Eq. (|3) by taking the correlation between 
filtered fields. For an isotropic survey covering the full sky with uniform noise, the minimum-variance lensing recon- 
struction has been worked out by Okamoto & Hu flF . We write their estimator here as 

^AlJ dnYl^,{h)V' [V{h)V,U{h)] , (4) 

where is a multipole-dependent normalization, and V{n) and C/(n) are the filtered maps, given as 

LM ^,expt 

= E Td^^i^YUi^)- (6) 

Here, are the multipoles of the measured sky temperature (after correcting for the instrumental beam) and 

C^^pj is the ensemble-averaged power spectrum for the measured temperature, nominally Cf^ + i where Cf^ 
is the (lensed) CMB power spectrum and Nf^ is the power spectrum of the experimental noise. In this work as a 
notational aid we will tend to use L to denote a multipole of the lensing potential, and t to denote a multipole of 
temperature. The power spectrum Cf^ in the filter for the J7(fi) field is of the unlensed temperature. This estimator 
has a nice interpretation. It is a correlator of the inverse-variance-weighted temperature map V^(n) with the gradient 
of the Wiener reconstruction U{i\) of the unlensed CMB, which agrees with the intuition above. The inverse-variance 
weighting reduces the magnitude of the Gaussian noise, i.e. chance 0V0 correlations, which contributes variance 
to the estimator. The divergence operation in the integrand of Eq. Q serves to extract the gradient part of the 
estimated deflection field. 

The approach of Eq. Q , has already been used successfully to find evidence for lensing effects in the WMAP data 
by cross-correlation with external data sets [171 HH] ■ The next step in the study of CMB lensing will be an entirely 
internal detection. It is expected that this will be possible in the near future with data from the Planck satellite. 
Planck noise levels should enable us to measure the lensing power spectrum with S/N > 1 over the multipole range 
^ =10-300. This will prove cosmologically interesting. For example, massive neutrinos should have percent-level 
effects on the lensing power spectrum if the sum of their mass eigenstates is ~0.1eV [19]. This is close to the lower 
limit on the summed neutrino masses inferred from oscillation data, and so Planck will at the very least be able to 
place useful upper limits [T5] . 

For the purpose of obtaining internal parameter constraints we require not (p, but an estimate of its power spectrum 
Cff . A straightforward approach to obtaining this is to form the estimator 

I ^ . . 

CV' = 2L + 1 E ^LM<i>*LM- (7) 

Apart from the well-known bias in this estimate associated with the statistical noise in the reconstruction of 0, there 
are several intrinsic complications which are usually neglected. The first is that although Eq. Q recovers the lensing 
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potential (pLM when averaged over CMB realizations, the actual reconstruction for a single realization is sensitive to 
the lensing potential over a variety of multipoles. In the power spectrum reconstruction of Eq. ([7| these modes may 
couple to give additional biases. The magnitude of these terms was first pointed out on the fiat sky in [5D]. 

A second and more general point of concern is that Eq. Q only reconstructs (I)lm to first order. The lensing 
expansion of Eq. ^ does not generally converge as fast as one would expect given the small size of the lensing 
deflections, and modern observations already test its accuracy. An accurate calculation of the lensed power spectrum 
Cf®, for example, was given in [15] . There, it was found that working only to leading order in Cf^^ [i.e. retaining 
the first two non-trivial terms in the Taylor expansion in Eq. ^ which is necessary for a consistent calculation of the 
power spectrum to 0{C'ff)] is only accurate at the 5-10% level relative to the true lensing corrections for multipoles 
£ =1000-2000. It is quite reasonable to find a partial breakdown of the low-order Taylor expansion in this regime 
— at these multipoles the typical angular scale of the CMB fluctuations is 5-10 arcmin, whereas the magnitude of 
the deflection |V(/)| is typically 3 arcmin. Things are not quite as bad as these arguments would suggest, since it 
is only the relative displacement of points that matters for the statistics of the lensed CMB, but it is still the case 
that there is general need for care when working with the lensing expansion of Eq. ([2]). Indeed, some numerical 
work has already been carried out [5T] which shows the effects of higher-order terms from Eq. ([2| in flat-sky power 
spectrum reconstruction, although they only found serious biases introduced at high multipoles where astrophysical 
non-Gaussianities also begin to gain importance. 

Finally, related to both of the first two issues, there exist intrinsic correlations between the multipoles of the 
reconstructed power spectrum, even on the full sky. These are to be expected since the 4> reconstruction is non- 
Gaussian, being quadratic in the lensed CMB temperature (which is itself mildly non-Gaussian). The magnitude of 
these correlations was investigated in the flat-sky limit in [20]. They found negligible correlations of 0[< 10~^) for 
all multipole pairs inside L < 2000. Their calculation treated the lensed CMB as Gaussian, however, and numerical 
verification that the non-Gaussianities introduced by lensing (which couple much more readily) are negligible would 
be useful. 

In this paper, we use numerical simulations backed up by analytical work to check that the current understanding 

of the reconstruction Cff' is adequate. We do this for the simplest possible case: on the full sky with homogeneous 
noise and no secondary effects other than lensing. Our simulations are discussed in Sec. [IT] In Sec. |III| we discuss power 
spectrum biases. Our principle new result is that higher-order terms from the expansion of lensing effects significantly 

suppress power at low-L in the reconstruction of C^"^. We provide simple analytic approximations for the bias which 
capture the majority of the effect seen in our simulations. In Sec. IV we give our results on the covariance of the 
reconstructed power spectrum, and derive further simple formulae which accurately describe the covariances on the 
full sky. Finally, we conclude with a discussion of the cosmological relevance of all of these complications, and discuss 
improvements to the usual lensing estimator which can largely remove them. Appendices provide technical details on 
the harmonic-space form of the lensing Taylor expansion and a flat-sky calculation of the main results in Sec. |III[ 



II. SIMULATIONS 

Our simulation procedure is as follows. We start with Gaussian realizations of (f> and 8 from spectra for a ACDM 
cosmology with /i = 0.73, = 0.022, 0^/1^ = 0.127, =0.95, (78 = 0.743 and three species of massless neutrinos, 
up to ^inp = 3000. Our lensing power spectrum Cff , for example, is shown in Fig. [l] We neglect the cf^'^ correlation 
induced on large scales by the integrated Sachs- Wolfe effect. The \ow-£ temperature multipoles play a negligible part in 
the lensing reconstruction and this choice has no effect on our results. We proceed to lens our 0(n) realization using the 
LensPix framework |22j . This takes the gradient of cf) using a fast spin-1 transform and remaps the temperature with the 
resulting deflection field, using cubic interpolation to approximate Eq. ([l]) accurately. We obtain a lensed HEALPix 
map at resolution Ngidc — 2048, which we convolve with a Gaussian beam of full- width at half- maximum ctfwhm — 
7 arcmin. We then add uncorrelated Gaussian pixel noise with standard deviation crjv = 27 /iK-arcmin. Finally, 
we deconvolve the experimental beam and harmonically transform our map back into multipole space, dropping all 
multipoles past ^^ax = 2750. We drop higher multipoles because the lensed power spectrum Cf ® results from a 
convolution in multipole space of width AL < 250, and so the lensed multipoles simulated with an i\n-p = 3000 cutoff 
are not cosmically accurate above this ^max- We have verified that for the noise level which we use the contribution of 
information in these high multipoles to the lensing reconstruction is sufficiently small that the location of our cutoff 
has negligible effect on our results. The power spectra of our lensed maps agree well with those expected, namely 



tCMB 



where Cf® is the lensed power spectrum theoretically predicted by CAMB 24 , using the prescription of [15) . 



4 



For lensing reconstruction, we perform Eq. (|4|. The gradient operations are again made with spin-1 transforms. 
The filtering operations are performed in harmonic space, with Cf^ our input unlensed power spectrum and Cf'^^^ 
as given above. 



III. POWER SPECTRUM BIASES 

To determine the expectation value of C^"^ it is most natural to work in harmonic space. Here the lensing expansion 
becomes somewhat more complicated, but the underlying idea is the same as in real space. Schematically, we have 

Qim = Qim + seem + S^Qirn + S^^lm + ' ' ' , (9) 

where the power of 6 denotes the order in (/>. Simplified expressions for the expansion terms are given in Appendix [A] 
In harmonic space, the Okamoto & Hu [16| estimator becomes 



where the symbol in braces is the 3j symbol and from here onward we use the compact notation that — m^}. 
The weight function ge^^ir^ (L) gives the filtered maps and performs the gradient operation: 

where the function F^-^li^ describes the rotationally- invariant part of the coupling between the three multipoles 
{ii,£2,L) and is given by 

Fe,Le2 = [El-El+El]Ue,LeAlQn)-'^'[^iQ^S) ■ (12) 
Note that Fi^^i^ is symmetric in its last two indices. For later convenience, we have introduced the compact notation 



-■a- ■■n 



^a{a + l)---n{n+l) (13) 
T^a-n = V(2a+l)---(2n + l). (14) 

In terms of these, the estimator normalization is 

A'l' ^ ^^fttU29hi2{L), (15) 

which ensures that the average of 4>lm over the CMB fluctuations /or fixed lenses is the true 4>lm- 

Averaged over realizations of the CMB and lenses, we can now see that our reconstructed power spectrum has an 
expectation value given by 

4li ^\M f h h L \ ( Ij, li L 



Thus, we need to determine the 4-point function of the lensed CMB. The dominant term is from the disconnected 
part, which is derived under the assumption that the lensed CMB is a Gaussian random field on the sphere: 



A 



L ~ E[5^i^2(^)]^[2Q^^xptQt^xpt] = ^L, (17) 



where the last equality holds only when the weights are chosen (optimally) as in Eq. (111. Here we have taken up 
the notation that ivj^''' gives the power spectrum bias with order p dependence on the lensing power spectrum Cff , 
exempting terms where it enters as a component of the lensed CMB power spectrum. 
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FIG. 1: Simulations of lensing (deflection) power spectrum recovery . The (+) symbols mark the average power estimate 
for 1000 simulations, and the gray boxed regions indicate the standard deviation for one realization (with binning over the 
multipoles corresponding to the width of the box). The black solid line is the input lensing spectrum and the dash ed lines are 
analytical calculations of the bias terms Nf' (black), N 



(red) and N^^ (orange). These are derived in Sec. Ill 



In Fig. [l] we show simulated reconstructions of the lensing power spectrum. We can see that the N^^ noise bias 

is the dominant source of error in the reconstruction. The noise bias in the deflection power spectrum L{L + l)Cff' 
is approximately scale-invariant for L < 100. We can understand this behaviour by noting that the dominant 
contribution to N^^"^ is from chance correlations between modes near the resolution limit of the experiment. For low 
L, we can therefore approximate L <^ £i, 
Taylor expand the C^^ power spectra in c 



2 in Eq. (17). Provided L is small compared to the acoustic scale, we can 



, {L) about £i to give 



->ee 

'2,expt 
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£i,cxpt 



r2 , 
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4E^ 



dim 



(18) 



Summing over £2, the leading-order terms for £1 3> i give 



1 .cxpt^j 



ee 

2,CXpt 



rr2 ■=■4 



^"-"^i^oxpt^ 



dlnCf^ 3 / dlnCf^^ 



dlnli 



d\n£i 



(19) 



Substituting into Eq. ([l7| and summing over £1 gives [L{L + l)]^ivj°''/(27r) = 2 x 10 ^ in good agreement with the 
large-scale plateau in Fig. [l] The rough scaling with sensitivity and resolution is 



[i(L+l)p<V(27r)^l/e, 



(20) 



where, here, ^max is an effective maximum multipole below which the measurement of the CMB power spectrum 
is cosmic-variance limited. Note that ^^ax defines the effective number of modes and that the noise power scales 
inversely with this number. 

While Nf^ gives the dominant bias, higher-order terms do have significant effects. This is illustrated more clearly 



in Fig. ^ 



where we plot C^' 



N 



(0) 



G1 . We therefore now set out to investigate the non-Gaussian terms which 
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arise from the connected part of the 4-point function — (OfjOf^Gf^Of^);^- To calculate these terms, it is useful to 
have a system for encoding the constraints which symmetry places on their form. The formalism for doing this has 
been worked out by [5S], and we briefly review the notation here. 

Statistical isotropy demands that the connected 4-point function function takes the form 
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LM 



( h i2 


-I. 




i4 




I nil 1712 




\ TO3 


7714 





(21) 



The term T^^^^ (L) is known as the trispectrum. Inserting this template for the 4-pt function into Eq. ( 16 1 leads to 
the simplified expression 



(22) 



ci ^2 14 



The trispectrum contains additional symmetries due to our ability to reorder £1 ■ ■ ■ £4. These may be explicitly encoded 
by taking 



i-iy 



£1 £2 L 
£4 £3 L' 



(23) 



There are eight remaining symmetries in the partially-reduced trispectrum representation P^^^^'^ (L) — reorderings 
within the pairings and exchange of upper and lower indices. We may again explicitly symmetrize with respect to 
pair reorderings via the reduced trispectrum T, such that 



(24) 



where Sc/ = £1-1-^2 + ^ and = £3+£i + L. Finally, we may satisfy the remaining upper/lower symmetry with the 
fully- reduced trispectrum T^^^^ such that 



T' 
If 



(25) 



The quantity T^^^^ (L) is an arbitrary function of its arguments, and its use enforces all of the symmetry requirements of 
the trispectrum. These considerations greatly simplify the computation of the lensed trispectrum from the expansion 
of Eq. ([9]), as they make it straightforward to differentiate between terms which are related by symmetries, and terms 
which are genuinely independent. 

In terms of the fully-reduced trispectrum, we can write 



(6,^6,^6,36,^0 -^E 



LM 





i2 




(£3 h 






TO2 


-M J 


I 7713 rn4 


M ) 



(-l)*^T^,^Jj(L)+ perms.. 



(26) 



which includes a sum over all 24 permutations of £1, £2, and I4. This form is particularly convenient for theoretical 
calculations of the trispectrum from given physical effects. It should be noted that fully-reduced trispectrum is only 
specified up to the equivalence class that gives the same trispectrum T^^^^(L) under the operations in Eqs. (23)-(25|. 
Expanding the trispectrum in Eq. ( 22 ) in terms of the fully-reduced trispectrum and using the symmetry of the 
weights 5,1,2 {L) and that they enforce e.g. £1 + £2 + L = even, we have 

'4A? 



E 9iii2{L)gi:,iAL) 



«1 1^2 t3 t4 



(27) 



With this considerable array of notation behind us, we are now in a position to calculate the effects of non-Gaussian 
couplings on the value of (C^"^). 
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A. 0{(f)'^ 

At second order in the connected trispectrum term is due entirely to correlations of the form {SQi^SQe^Qi^QeJ . 
It has the fully-reduced form [25] 

^^1^2/7^^ _ f^4>4> f 



Inserting the primary {Pi^g^ ) form of this term into Eq. ( 22 ) (or using Eq. ( 27 1 directly) results in 



{CV')c = I^Y9iML)fi^Li2 ] (primary term). (29) 



For any choice of weights, this reduces to which it was our intention to reconstruct, provided that the estimator 



is normalized so that the map-level reconstruction is unbiased, i.e. Eq. (15 1 is satisfied. The normalization requires 



knowledge of the unlensed CMB power spectrum and any error in the fiducial spectrum used leads to a multiplicative 



bias in the reconstructed power. Inserting the other pairings into Eq. (22 1, on the other hand, results in the first-order 
power spectrum bias N'j^^: 

ivi^'=2('^y Y i-lY''^''9i^i2iL)gi,iAL)lllf^^r\ctfh,L'i,^ (30) 

^1^2^3^4-^^' 

This N^"^ term is the all-sky analogue of the Nqq ee(-^) noise bias of Kesden, Cooray, and Kamionkowski [20] (KCK 
hereafter; see also Appendix [b| . In practice, we find that the flat-sky result gives a good description of the biases 
found in our simulations. This is a reasonable result as the magnitude of the N^^ bias is non-negligible only at large 
multipoles and here flat-sky calculations converge with their harmonic counterparts. For the remainder of this paper, 
we will calculate N^"^ using the flat-sky expression, with all-sky power spectra as input. 

In Fig. [l] we plot the N^"^ bias. At low multipoles it is an insignificant contribution, but at L>250 it produces an 
0(10%) effect. At L> 1000 it begins to dominate over the Cff' term which it was our intention to reconstruct. We 

show in Fig. [2] that at high multipoles the N^'^ bias provides a good fit to the excess power seen in our simulation 
results. 

There are a number of ways to account for this behaviour. Hu showed that for lensing, the primary {Pf^f^) 
component of the trispectrum dominates over the other couplings when min(£i, i!2j ^3, ^4) 3> L [25j The primary terms 
represent flattened quadrilaterals with the short diagonal of length L supported by lensing, the power of which is 
strongly peaked at low multipoles. A heuristic (but frequently useful) way of rephrasing this argument in a way 
which is speciflc to lens reconstruction comes from considering the effect of the weights gg^i^{L). For ^1 3> L, the 
weights are only non-zero for £2 in a narrow range of width 2L centered on ^l (due to the triangle constraint). The 
term imposes couplings between two weight terms, and so the region of (^1,^2,^37^4) space over which it can 
accumulate signal is correspondingly suppressed. In general, at low-L the dominant trispectrum terms which bias 

Cff are therefore those which factor maximally under the weights. At high-L, which terms will dominate is less clear. 

The eventual dominance of the N^'j^^ bias was explained by KCK as the interference of lensing modes at L' < L on 
the reconstruction. It appears that in this case the cross-coupling nature of the non-primary terms which suppresses 
their contribution at low-L allows them to dominate later on. As C'j^'^ is very red, the coupling of the N^'^ term to 
C^*^ power at lower-L allows it to increase relative to the primary term. 

B. 0(9!)*) terms 

At second order in the power spectrum, there are several fundamental groupings of 5"Q and unlensed terms in the 
lensing expansion which contribute to the fully-reduced trispectrum. These are all contained in four fully-reduced 
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trispectra: 



5e(5e<5e<5e 



4 4 L 

4 4 4 



E 



16^ 



f-<BQ r^BQ r^4>4> r<'t 



rr2 ■='2 



, r"f"f' T? ST^ r<^^ r"f"f' ( p \ 



ii h L 

£a 



s^eseee 



("£1 - 2/3) Cf^^Fi^Lti 
Here, the value R is half of the mean-square deflection 



R 



E 



An 



(31) 



(32) 



and the (g)^^^J(L) function is given in Appendix A The summation Yl'i^e,, the second grouping is restricted to £a and 
£b such that 4+^6+^1+^2 — even and ^a+^b+^3+4 = even, thus enforcing the parity constraint ^1 +4+^3+4 — even. 
Similarly, the summation in the third grouping is restricted such that £3 + + ia + it — even. 

Many of these terms are difficult to calculate numerically. Fortunately, the dominant subset is calculable: 



dom 



F, 



I2LI1 



Y^ct^Cf^iF,,, 



■=■2 I T'2 T'2 
rr2 ■=■2 



2R{Zi - f /3)Q« 



ee 



(33) 



The choice of the terms in Eq. (33 1 is motivated by the discussion of the previous subsection, in that these terms are 



the subset of Eq. (|31j) which maximally factor under the weights for the primary trispectrum coupling. Calculating 

dom 

the effect of T^^^^(L) through the primary trispectrum coupling results in a \ow-£ bias at 0(0'*): 



■=-2 I T'2 ■=■2 



rr2 ■=■2 



2i?(El - 1/3)C« 



ee 



(34) 



We plot these dominant terms to the n'^j^^ bias in Figs. [1] and [2] Overall, they contribute a negative bias to the 
reconstructed power spectrum which is w 15% of C^*^ for the concordance cosmology which we use. They provide a 
good fit to the negative bias seen in our simulations for L < 200. 

The dominant terms that we have retained in Eq. ( 34 1 share a common origin: they are due to correlations of 
the gradient-order potential reconstruction of (j)LM with higher order terms. Equivalent terms arise, for example, if 




FIG. 2: Biases in . The points are the average of 1000 simulations and the boxes give the measured error in the recon- 



struction for one realization (with binning over the multipoles corresponding to the width of the box 



Sec. 



Ill 



are also shown: N'^''' dotted (red); A*'^^' dashed (orange); and their sum in solid (blue). 



a good fit to the observed power in simulations. 



The biases calculated in 
These calculated biases provide 



one calculates {Cff') — the cross-spectrum between the input and reconstructed lensing potentials. If we write the 
reconstructed (f>LAi as^ 

(t>LM = n% + <PLM + n^LM + '^LM + ^fll + ' ' ' ' (35) 

where, for example, tS^^j^j is the statistical noise in the reconstruction that is second order in and arises from Q5^Q 
and 5Q5Q terms, then the 0{Cff') contribution to {Cf ) is from 4'lm only (since (n^^2^)cMB = by construction) 
and the 0[{Cff')^] contribution is from nflj. Correct to fourth order in 0, we find 

dom 



(Cf)=Cf + -ivP, (36) 



dom 

where A^^ — '^{4^):,m'^lm) given by the right-hand side of Eq. (|34|) with one factor of the normahzation ehminated 
with Eq. ( [151 ). 

(2) 

We can develop a useful approximation to N)^ on large scales by noting that the sum over ii and I2 is dominated 
by small-scale modes near the resolution limit of the experiment. Since C^"^ is so red, £1, ^ £a ^ f-i in the remaining 
summation and we can approximate 

W„ i^M^eJ =^2712 -^Z^^^a — fj2 ■ (.37) 



^ We are grateful to Antony Lewis for suggesting this approach. 
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Noting further that the second-order result for the change in the temperature power spectrum due to lensing is 



(38) 



we can approximate 



N 



(2) 



4C 



00 



ni 



ee 



(39) 



(2) 

The fractional Nj^ bias is therefore largely determined by the difference between the lensed and unlensed power 



spectra. We can use similar methods to those which gave Eq. ( 19 1 to approximate Nj^ further at low L 



N 



(2) 



47r 



n 



cP®Ac; 



ee 



1 



2dln^ 



ee^j 



3 dCff dACf^^ 
8dln£i dln^i 



(40) 



where we have used Eq. ( 15 1 for the normalization. The only L dependence is in the prefactor, but Eqs. ( 17 ) and ( 19 ) 

r(2) 



show that A^'. 



const, at low L so the iV, 



bias is a constant fractional bias there. Numerically, the dominant 
contribution to Eq. (|40|) is from the final term and the bias evaluates to —20%, in reasonable agreement with the full 



evaluation of Eq. (34) reported in Fig. [2] The bias become less important for surveys with lower resolution: for an 
experiment limited to a maximum multipole ^max = 750, we find the bias reduces to —4%. 

We have seen that the dominant terms in the 0{(f>'^) trispectrum provide a good explanation of the low-£ bias seen 
in our simulation results. Figure [2] provides empirical evidence that h igh er-order terms do not bias the reconstructed 
power further by more than ^ 3%. This seems plausible given Eq. (39). Replacing AC®® by its non-perturbative 
equivalent (see e.g. Ref. llSj), the change in AC®® is around 20% at 



£i ~ 2000. To the extent that this replacement 
mimics the non-perturbative generalization of the low-L bias (see Appendix [b] for further justification), we should 
expect higher-order terms to produce a fractional bias ^ 0.15 x 0.2 = 0.03. 



IV. POWER SPECTRUM COVARIANCE 

We now proceed to investigate the covariance of our power spectrum estimates. A similar calculation has been 
performed on the flat sky by [5D], however the full-sky calculation is somewhat cleaner. 
We are interested in the covariance of the reconstructed power spectrum which is given by 

Cov(C^'^,C^,'^) = {cf'Cff) ~ {Cf^^){Cp). (41) 



The cross term. 



X 



ctctf}-i^ E E E E i-^y 

li l2 L\f U L L'\f l^ 4 V 

mi 7712 M J I 7773 ^'i — A-f J I 7775 7776 M' J I 7777 ^778 —M' 

X 5,^,, {L)g,,,, {L)g,,,, {L')g,,,, (L') (6,^ 8,^ 6,3 8,^ 8,^ 8,^ 8,^ 8,^ ) , (42) 

requires us to calculate the 8-point lensed correlation function. Since the lensed CMB has zero mean, and we are 

ignoring the low-^ Cf'^ correlation and hence all odd connected correlations, the 8-point function can be expressed 
in terms of the connected 2-, 4-, 6- and 8-point functions. The latter two of these are at least O(0^) 27^ and so we 
neglect them on the grounds that they should be sub-dominant. 

We first consider those contributions to the 8-point function that involve the trispectrum. There are 70 terms 
involving products of two trispectra and a further 210 involving the product of a trispectrum with two 2-point 
functions. Guided by our experience with the power spectrum, it seems likely that the dominant couplings are those 
which factor most under the weights. This motivates us to consider only those terms with trispectra which couple 
terms between two weights, e.g. {Qi^Qi^Qi^Qt^c ■ There are three such terms involving a product of trispectra and 
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18 involving a single trispectrum. It is convenient to re-express these in terms of the full 4-point functions and their 
Gaussian (disconnected) parts, so their contribution to the 8-point function is 

{QeeiOeQiJiOieiee^ei,) - ((e,^e,^)(e,,e|j + 2 perms) 

X + 2 perms) + (34 o 56) + (34 o 78). (43) 



The 4-point terms make simple contributions to Eq. (42): 



(e,^ e,^ Qt, Qt, ) {Qt, Qi, Qi, Oi, 

For L = L' this leads to a dominant variance of 

Var(C^^) = 



5LL'{Ct^){ctf)/i2L+l) 
5LL'{Ct^){ctf)/{2L+l). 



2L 



(44) 



(45) 



which agrees with the results of [20] on the flat sky. It also justifies the assumption of |28j that the statistical noise in 
the (j) reconstruction behaves predominantly like Gaussian noise giving the usual cosmic-variance result for the power 
spectrum variance. This is a slightly different version of cosmic variance than we usually deal with, however, in that 

the terms in C^"^ arise from correlations between a large number of multipoles. 



We now consider the disconnected part of the 8-point function and the remaining 2-point terms in Eq. ( 43 1 . As we 



shall see, these generate a covariance between the power spectrum estimates. The disconnected part of the 8-point 
function involves 105 terms but 45 of these are zero (for L and L' ^ 0) because they correlate one or more pairs of 
Q which are jointly weighted, which isolates {4>lm) = when averaged over realizations of large scale structure. The 
remaining 60 non-zero terms are characterized by four fundamental pairings, with terms within each pairing giving 
equal contributions to the covariance by symmetry: 



e. e, e, e, 



e, e, e, e, 



e. e. e, 



e. e. 



9,,9,^9,^9,^ 



9. 9. 9. 9, 



9. 9, 9. 9. 



(4 terms) 

(8 terms) 
(32 terms) 

(16 terms). 



(46a) 

(46b) 
(46c) 

(46d) 



The first and second pairings cancel with the 12 non-zero 2-point terms in Eq. (|43j). The third pairing results in a 
covariance of 



{Ct^Ctf)a = 32 



The final pairing is more tightly coupled, and results in 



L I2 



,GXpt 



expt 



(47) 
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n 



LL' 



L+L' 



L (-1 I2 \ ^ee ^ee ^ee ^ee 

L' £3 £4 I ^l'CXpt'-'€2,OXpt'-'f3,OXpt^^4,CXpt 



X 5,,,, {L)ge,e, {L)gi^i, {L')gi,t, (i') • (48) 



We expect that the terms of Eq. ( |47[ ) should be dominant as they factor most under the weights, and so from here on 
we consider them only. In Fig. [3] we plot the correlation matrix 



n{L,L') = 



Cov(cf^,Cy) 
Var(C^^)Var(C^,'^) 



(49) 
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The agreem ent between the measured and theoretically-approximated correlations is excellent. Using the approxima- 
tion in Eq. ( 19 ), the correlation matrix for low L and L' should scale as 7?.(L, L') ^ %/ LL' H\-^^y^ in agreement with the 
Note that the correlations decrease as the resolution of the experiment increases. 



arguments in Ref. EOj. Note that the correlations decrease as the resolution of the experiment increases. While the 
shape of our numerical covariance agrees with |20j , the magnitude we find is at least one order larger. The correlations 
are at a level of < 1%, although we note that this correlation is for unbinned spectra. With such broad correlations, 
binning increases the correlation roughly in proportion to the bin width. The binned spectra in Fig. [T] for example, 
have correlations of 0(1%) for i < 1000, where logarithmic binning is used, and of 0(10%) for i > 1000, where the 
bins are rather wider. In the discussion which follows, however, we will show how these correlations may be almost 
entirely eliminated in practice. 



V. DISCUSSION 



We have identified three complications with the usual lensing estimator: 

1. the 0((/)^) n'i^^ excess power bias; 

2. the 0((/)'^) A^l^"* power suppression; and 

3. a small intrinsic covariance between the estimated Cff' . 

We will discuss the implications of each of these terms and their remediation in turn. 

We begin by investigating the relevance of the bias terms for cosmological parameter determination with a Fisher 
matrix approach. Our implementation follows closely that of and we refer the interested reader there for details. 
The incorporation of C^"^ reconstruction into a parameter analysis allows one to produce dramatically improved 
constraints on the 'dark' parameters which affect the late-time evolution of the Universe [1^1 . With Planck for example, 
a factor of two improvement over what is achievable without lensing reconstruction is forecasted for the determination 
of the (summed) neutrino mass in simple models, and rather more in models with additional parameters such as 
dynamical dark energy [T31 [21] • We begin by asking which region of the Cff' power spectrum this determination is 
made from. We use the fiducial cosmology given in Sec. plj but take fl^h'^ = 0.006 (i.e. X^i/™!^ = 0.6eV). This is 
large enough that it would be detected with significance hyPlanck , which is important to prevent the hard n^h^ > 
prior from corrupting the Fisher analysis. In the top panel of Fig. [i] we plot the relative constraints on il^h^ for our 
simplified version of Planck , assuming that the reconstructed Cff' power spectrum is ignored above some value of L. 
We can see from this figure that for Planck , most of the constraining power of lensing on neutrino masses comes from 
the multipole range 100<i<700, which will receive large contributions from both bias terms. 

We now consider (pessimistically) the effect of the N^^'' and n'^'' biases if completely uncorrected. We can relate 
biases in the 'observed' lensing power spectrum to parameter biases by perturbing the Fisher-approximated likelihood. 
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FIG. 4: Relationship between reconstruction and parameter constraints for Planck . Top panel: constrains on Q^h^ 
assuming that the Cff' reconstruction is discarded for L > £, divided by the constraint if the full reconstruction were included. 
Inclusion of lens reconstruction improves constraints on by a factor of two, consistent with e.g. [T3]. Bottom panel: bias 

in the determination of f2„/i'^, as a fraction of the random error, induced by the A^^^' (black solid) and N^^ (red dashed) biases, 
assuming that the biases have been completely removed for L > £. 



finding that 

A0, = Fr:' E [dC^"" /d0A[AC^''][Cov-']AB,CD, (50) 

AB,CD 

where 9i indexes a cosmological parameter, Fij is the Fisher matrix, AC^^ is the power spectrum bias and [Gov] is 
the covariance matrix of the measured spectra. The labels A, B, C and D run over the observable fields T, E and 
(f). In the lower panel of Fig. H we assume that the bias is uncorrected for L < £ and perfectly removed for L > £. 

(2) 

We can see that neglect of either bias leads to la errors in determination of the neutrino energy density. The iV^ 
bias suppresses the reconstructed power, particularly at low-^, and this partly mimics the effect of massive neutrinos, 
leading to an overestimate of fl^h^. This is particularly worrisome if we are placing an upper limit on the neutrino 
mass, as neglect of the Nj^ ' bias would lead to a spurious (albeit marginal) detection of massive neutrinos. We 
conclude that an accurate treatment of both biases must be made in order to obtain accurate parameter constraints, 
and so now turn to ways of mitigating these. 

For the N^j^^ bias, the suggestion for removal has been to perform iterative estimation and subtraction [201 

An additional approach, which is perhaps conceptually more straightforward, is to view N^"^ as a normalization 
effect which results in an estimator which is a convolution over the true lensing power spectrum rather than a direct 
estimate. The kernel of this convolution can be calculated and inverted to produce an unbiased estimate of Cff . 
We note that this kernel depends on the true unlensed CMB power spectrum. In practice, our uncertainty in this 
quantity therefore determines the extent to which the n'^j^^ bias may be treated. Due to the large noise with which 
upcoming experiments will reconstruct the lensing potential, both of these procedures are likely to be unstable unless 
we can take the shape of the lensing power to be characterized by some small set of numbers (such as the cosmological 
parameters of interest). 

Experimental realism introduces additional difficulties, however. In particular, we note that for a realistic ex- 
periment with inhomogeneous sky coverage, the normalization of a quadratic lensing estimator must generally be 
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determined with simulations [T7]. To treat the N^^^ bias using off-diagonal entries in a normalization matrix would 

require that we be able to Monte-Carlo these elements — a difhcult task. In this case, the N^'^ bias will most likely 
need to be absorbed into the normalization. Suppose, for example, that the normalization is determined with lensed 
realizations of the CMB using some fiducial Cff' spectrum, as in [17]. At each L, a normalization Al is determined 
from the ratio of the Monte-Carlo average of the reconstructed power (after correction for the A^^"'' noise) to the 
fiducial power there. Assuming that the shape of the fiducial C^*^ matches the underlying lensing power, this results 
in an unbiased estimator at the power spectrum level, albeit with slightly greater variance than would exist in the 
absence of the N^'' bias. The quality of our assumptions about the underlying C^"^ shape would then need to be 
quantified and folded into the systematic-error estimate. 

Treatment of the N^"^ bias presents further difficulties due to its non-linearity in Cff' . It is particularly troublesome 
in the case where the normalization is determined by Monte-Carlo simulations, as care must be made to determine 
the normalization over a small range of Cff" power, to ensure that the estimator response is approximately linear. 

(2) 

The Nl^' bias would then appear as an offset which could be subtracted in a fiducial model with simulations using the 

fiducial Cff" power. This would clearly be a cumbersome procedure. We find, however, that a small modification of 
the optimal estimator can effectivel y re move this bias prior to normalization. Our starting point is the approximate 



(2) 

expression for the bias in Eq. (37 1. In terms of this, we have 



c 



L I Hi 



c 



L 1^2 



\i1i2 / 



+ A^f +A^P 



1U2 



llt2 



+ Aff +A^P, 



(51) 



where /^^^^ = AC^^/^jL^^ -F ACf^^ fe^Lii and, recall, ACf^ is the change in the CMB power spectrum due to 
lensing. We now see that if we modify the weights to 



9iii2{L) 



(52) 



where f^^j^i^ is computed for some fiducial Cff" , the quadratic part of the estimator response to power Cff' becomes 
linear times the difference between the true lensed CMB power and that in the fiducial model. The estimator normal- 



ization, for an ideal survey, is still given by Eq. (17) using the original weights. In practice, the normalization could 



be obtained by Monte-Carlo with simulations using the fiducial Cj/ and the modified weights in the reconstruction. 
The cost of this modification is an increase in the estimator variance of ~ 15% for L < 300, however this is most 
likely justified given the improved bias properties of the estimator. Additionally, this technique makes clear that our 
ability to debias the estimator is determined by our understanding of Cff' only through the lensed and unlensed CMB 
temperature power spectra, a simplicity which one would not necessarily expect for a higher-order bias. This insight 
makes the assessment of systematic errors due to uncertainties in cosmology much more straightforward. 

Another, more straightforward method of treating the bias is motivated by the recent work of Lewis et. al. j30| . 
who point out that to a good degree of approximation, the response of the lensed CMB covariance to changes in the 
lensing potential is determined by lensed power spectra. This is in contrast to the first-order formula used to derive 
the quadratic estimators of Okamoto and Hu |16| . which contains an unlensed spectrum. Incorporating this insight 
into the quadratic derivation one would obtain the same estimator, but with lensed rather than unlensed spectra 



in the filtering of Eq. (11). Because these appear twice in the estimator normalization (Eq. (15)), this is in fact 
equivalent at 0{(j)'^) to the correction described above. This method of debiasing is even simpler to implement than 
the one proposed above. 

Finally, we consider the covariance of the lensing power reconstruction. The loss of information which this covariance 
represents increases the effective error bars of the reconstruction, although we have not quantified this here as we 
find that an improved, internally-calibrated estimator does a good job of reducing this covariance. Similar ideas were 
put forward by Dvorkin and Smith in the context of optical-depth reconstruction [31 . For lens reconstruction, the 
modified estimator for the power spectrum is as follows. Instead of removing the N^^^ bias by direct subtraction, 



we subtract 2N'f^ — N^"', where NY'' involves the observed total power spectrum in our realization of the universe. 



r(0) 



r(0) 



15 



/vee . 



By construction, N'j^^ is quadratic in the observed temperature and has (N^/^^) = N'j^\ Subtracting 27Vj°'' — N'j^^ 
is very successful in reducing the intrinsic covariance of the estimator, while preserving the expectation value of the 

power. This is because the dominant terms of Cov{Cf^ , N^^}) + Cov(C^f , iv|°-*) are identical to those of Eq. (47 1, 

and the dominant terms of Cov(7Vj°\ TVj?'') are 1/4 of those in Eq. (47). Their combination therefore cancels the 



dominant contribution to the off-diagonal covariance from (C^ )a. With this improved estimator, we find that 
the covariance is reduced to a level which is unmeasurable in our simulations (< 0.05%). For an experiment with 
only partial sky coverage, N^"^ can be determined by simulations for each of which the quadratic (j) reconstruction is 
performed by correlating independent Gaussian CMB maps drawn from Cf^^^ and Cf^^^ respectively. 

As we move toward experimental realism, the number of quantities which are determined by Monte-Carlo becomes 
increasingly worrisome, although some sanity checks are available. The accuracy of the n'^^ subtraction, for example, 
may be assessed by a jackknife test in which one reconstructs the lensing power using two estimates: one which 
correlates only even multipoles, and one which correlates only odd niultipoles. For parity-symmetric sky coverage, 
this results in an estimate of Cff which is free from the n'^-' bias, providing a useful consistency check. 



VI. CONCLUSIONS 



We have thoroughly investigated the behaviour of the optimal quadratic lensing estimator, both analytically and 
with simulations. We have discovered a new bias in the reconstructed power spectrum, and presented the lensed 
CMB trispectrum at 0{4>'^) to explain it. This bias has a physical interpretation as an anti-correlation between the 
first-order reconstruction of ((> and higher-order terms, and so it will also be relevant for cross-correlation studies. We 
suggest that the following estimator provides a good basis for the reconstruction of Cff': 

= Ai' Cif -2Nf ^ + N^°^ , (54) 

where Off denotes the usual optimal estimator of Eq. ^ with a slight modification of the weights, given by Eq. ( [52| . 
This estimator has a slightly larger variance than the standard optimal estimator at low L, however it is effectively 
free from contamination by higher-order lensing terms and has negligible covariance. Uncertainty in the underlying 
cosmology leads to potentially imperfect removal of the biasing terms, however we find analytically that this possibility 
is completely characterized by our understanding of the lensed and unlensed temperature spectra. 

In the near future, CMB lensing will make the transition from detection to precision science, and the concerns 
which we have addressed will be increasingly important. The reconstruction of the large-scale lensing potential is a 
demanding task, however the new window which it will open onto the contents of the Universe will most certainly 
provide a worthwhile view. 
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^ This modification to the estimator is also what one finds if one constructs a quartic estimator in the observed temperature for cf 
based on an approximate maximization of the likelihood for the lensed temperature (truncated at trispectrum order). 
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Appendix A: Simplification of Expansion Terms 

We wish to calculate the connected part of the trispcctrum to C'((/)"'), which requires a third-order expansion of the 
lensing effect. In harmonic space, these terms are given by 

SO em = J2i e 4'ii^e2^u^l"'^ 

s'eim = '^tli^eJi^hA.JriT"'' (Ai) 



s:3r\ 1 J. J. J, r\ T^mmim2m3m4 



where 



We may simplify the covariant derivatives in Eq. ( A2 ) by exploiting their relationship with the spin-raising and 
lowering operators [33 . Their application thus results in the generation of spin-s spherical harmonics, which are 
defined by the application of the spin-raising and lowering operators to the ordinary (spin-0) harmonics |34l I35j . 



Following the derivation in 16J, it is straightforward to show that 

(^r,"eV"+)e("-)) ^-ie^)eL"-)(V(^-.)(^ + . + l) [.+1^,™] el - ^(^ + + 1) l_,Yr] e^), (A3) 

where n_ — n+ = s, the (null) spin basis is given in terms of the unit polar and azimuthal basis vectors, e± = {eg±ie^), 

and we have defined e^-* = e^±e± ■ ■ ■ e^± ■ Repeated 
proceed to simplify the expansion terms individually. 



and we have defined e^-* = e^±e± ■ ■ ■ e^± ■ Repeated applications of V reduce to applications of Eq. (I A3|) . We now 



^ Tmmim2 

The simplification of Ilg^^"^^ has been given in the general case by |25| . Expanding the gradients of the spherical 
harmonics and then contracting indices gives 

-'«i«2 -'^2"fii2 a'^iye [1^11 -i^£2 + -i^fi i^<?2 )■ 

This vanishes by parity unless £ + ii + £2 — even. Now we may use an identity for integration over three spherical 
harmonics |36j . 

dn^Yr sjiT sjr:' = i-ir+^^ {[ ^i ] (± i\ ii) , (as) 



/47r V ~'*2 J \—iTi mi m2 
and the following result from the recursion relations between the 3j symbols |36j : 



£1 £2 4^ _ 1 f£i £2 £ 
1 -1 



'='2 T'2 ■Z'2 



to show that 



Iu,i2 -(-1) [-mm, 7712)^"^'^' 
where the parity constraint is enforced by F^^i^. This simple result for IJ^^^"^^ also follows rather more directly by 



repeatedly integrating by parts in Eq. ( A2 ) 



It follows from parity that J^^^^^"^^ vanishes unless £ + £, + £2 + £3 = even. When this is satisfied, we have 

^«Xr"' = =^V(4 + 2)(^3-i) / df7rr*_iy™^_ir,72r,7 + ^^i^|^ / dnvr^Y,^^ _,y,^\y,^\ (A8) 
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Thus the gradient integral reduces to an integral over four spin spherical harmonics. We may reduce this to a 
three-harmonic integral by invoking the Clebsch-Gordan expansion for spin harmonics ^37] 



Application of the three-harmonic integral identity of Eq. ( |A5| then gives 



71^ 2_.^"^ -S2 Sj [mi 7712 -Mj 



JO V™i V™2 v™* — 

""si-^^l S2^«2 S3-^l3 S4^£4 ~ 



47r 
X (-1) 

This may be used finally to give two compact expressions for 
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1 + i'2 + ^3 = even: 
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(A9) 



(AlO) 



(All) 



where 



£ h L\ (£2 £3 L\„ 

-1 1 1 0-1 

4 L\ ( £i £2 L\ „ 



n 1 -1 



2-2 



?i ^2 L 
-I -1 2 



^(^3-1)^3 + 2) 



(A12) 



Note that G^^^^ (i) and (L) correspond to the two distinct ways of coupling the arguments of J (it is symmetric on 
the second and third arguments). They are therefore related by a 6j recoupling coefficient (for ^ + +^2 +^3 ~ even): 



(A13) 



In our calculations in the main text, we only encounter terms with (£ + £i+L) and (^2+^3+i) even. In this case 
we can simplify further by using the additional recursion result 



ili2 ^3 \ 1 i2 is 

11-2 2 I 



■=■2 ("='2 _ ■='2 I T'2 \ I '='2 (-^2 I T'2 _ T'2 

"^1 '^"^3 "^1 '^^^2) ^ "£2 '-"<;3 ^ "^1 "<'2 

Hw3v/(^3-l)(4 + 2) 



(AM) 



which holds only for (^1-1-^2+^3) even. Using Eqs. (A6) and (A14| we find two simple equivalent forms for Qf,^\_^{L): 



e53 (^) = 



{^L"Ly^eeie2e3"eit 



FLi2i3("\ - + "i) 



(A15) 



O TV^mmi m2rn,3 7714 



Using the tools introduced thus far we may produce a general expression for the term. In our 

calculations, however, we only encounter terms in which K is contracted on its second and third or third and fourth 
sets of indices. The first of these is 



E/ -i-imi Tj-mmi—mimsmi 



dVLYP 



(AI6) 
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Invoking the result that 

^V.yrV.y/"* = ^-^^9^,, (A17) 

m 

where gij is the metric on the unit sphere, and using 

y2yfcym ^ v'=(v2 + i)y/" = (1 - Ej) v'^y^, (A18) 

we find that 

E/ I \mi Tv-mmi— mim3m4 / -i -3^2 \ K^^i^^^i) -rmm^ra^^ fAlQ") 



A similar calculation for the contraction of K on its third and fourth indices gives 



V J Mi£2^2^4 877- Ulf-i ■ V ' 



Appendix B: Flat-Sky Trispectrum 



For small patches of the sky, the curvature of the sphere is negligible and we can represent the CMB in terms of 
Fourier modes rather than spherical harmonics. Geometrical terms in this flat-sky basis are represented as scalar- 
products, which frequently makes them easier to work with than the corresponding full-sky expressions. Here we 
present expressions for the flat-sky trispectrum following the Fourier conventions of [3]. 

The flat-sky trispectrum, r(li, I2, 13, 14), is defined by 



(e(ii)e(i2)e(i3)e(i4))c = {2ti)-H{\i + 12 + 13 + i4)T(ii, I2, 13, 14). 



(Bl) 



The trispectrum is permutation symmetric in all its vector arguments and rotational and parity invariance mean that 
it depends on five scalar parameters of the quadrilateral formed from its arguments. As with the spherical trispectrum, 
it is frequently convenient to work with a fully-reduced trispectrum, T 

1 



(e(ii)e(i2)e(i3)e(i4))c = 



(2^ 



2 5(ll+l2+L),5(l3+l4-L)T 



l^^^^j(L), such that 

(i) + perms. 



(B2) 



where we include the sum over all permutations of li, I2, I3 and I4. The fully-reduced trispectrum is then an arbitrary 
function of its arguments but any permutation of £1, ^2, ^3 and i?4 within the (12) (34) pairing generates an equivalent 
trispectrum. The flat-sky fully-reduced trispectrum is related to its spherical equivalent by |25j 



^esiS^) - -^^^L^^iii2i^U I 



4 4 L 




{L). 



At second order in 



the flat-sky trispectrum is 
T(li, I2, 13, 14) = 2^^'^'-^^^^|ii+i2| 



,l(li+l2)-l2 (I3 + I4) ■ I4 + perms. 



Breaking up the delta-function in Eq. (|B1|), we have 

(e(ii)e(i2)e(i3)e(i4))c = 



(27r)2 



5{h + h+ L)(5(l3 + I4 



L)Qf Qfcf (-L.l2)L.l4 



perms. 



from which we can read off the fully-reduced trispectrum: 



T 



1 



£i+ii) 



(B3) 



(B4) 



(B5) 



(B6) 



Making use of Eq. (B3), in the limit of large arguments [so that e.g. i'^ ~ £{£ + 1)], we recover the spherical result in 
Eq. ((28l). 



19 



At 0((f>'^), there are four types of term that contribute to the trispectrum: 



(5e(5ei5e<5e 
12,13,14) 



T(li, 12,13,14) 

nil, 12,13,14) 



T(li, 12,13,14) 



3 / [1 • (li - 1) 1 • (I2 + l)(li + I3 - 1) • (ll - 1) (ll + I3 - 1) • (I2 + 1) 

X C^i\ I C'lttia I ^ ^|?i+i3-i| + perms. 

i^ee^ee j ^^i^ . ^ . 1 1, . (i _ 1^ _ 13)1, . (i _ 1^ ^ 13)^^^,^ ,^_,^| + perms. 

ll • (ll + i3)cr I4 • (I2 + 14) 

+ qt+m^^r 14 • (I2 + 14) / (^1 • (ll + 13) [1 • (1 - 13)]' crcp,3i 

- iCZ"" I j^h ■ (1 + 13)1 • (1 + l3)li • (1 - 12)1 • (1 - l2)Cf .^|C|tt,3| + P™- 
ll ■ (ll + h)CZ^ I4 • (I2 + I4) + perms., 



where 



R = 



An 



L-'CrdL 



(B7) 



is half the variance of the deflection field. 

The subset of trispectrum terms that dominate the power spectrum of the reconstructed lensing field are 



dom 



T (ll, I2, 13, 14) = -c^^+uiRelcZ"" ll • (ll + i2)Q«^ I3 • (I3 + 14) 



+ I3 • (I3 + 14) I (^1 • (ll + I2) [1 • (1 - li)]^ Cf 



+ perms., 



(B8) 



where we have chosen this particular permutation to display for later convenience. The first term (plus its permuta- 
tions) generates a fully-reduced trispectrum 



T 



(^1^2) 
3^4) 



2^- 



(B9) 



Using Eq. ( B3 1 , this gives a full-sky trispectrum which for large multipoles reduces simply to the second term in 
id term on 



Eq. (33). The second term on the right of Eq. (B8I gives a fully-reduced reduced trispectrum 



^^;^la • MK ■ (la - ll)]^C£°Cn_^_,^| 



-2Ct*CZ% • Lli • L 



d^la la ■ ll 

(2 



(BIO) 



where I3 -|- I4 = L = — (li -I- I2). To write this in a form that manifestly depends only on ii, £2, ^3, ^4 and L for 
comparison to the full-sky result, we write C^i'^_^^^ = J d'^lbS{la + h — h)Cf^ and use the large-f expansion of the 
delta function (e.g. [IS]), 



4.U f^a ib il 







ia ib il 

rua mb -mi 



(Bll) 



where, for example, is the angle that li makes with the x-axis. Substituting into Eq. (BIO), integrating over (pi^ 
and 0^ and simplifying with recursion relations for the 3j symbols, we find 



L')iil~ij^L')Y^CZ^C, 



eer«t><t> (^1 + ~ ^fc) 



(B12) 



This generates a full-sky trispectrum which for large multipoles reduces to the first term in Eq. ( 33 1 . 
For lens reconstruction, the standard optimal estimator is given on the flat-sky by |28j 



(L) = Al 



e(ii)e(i2)ffi,i,(L), 



(B13) 



20 



where L = li + I2 and 



5iii2(L) 



/hi.(L) 



^'-'£i,cxpt^£2,oxpt 



/i,i,(L) = li.LQ^e + l2-LC7| 



The normahzation is 



/iii2(L)giii2(L) 



(B14) 
(B15) 

(B16) 



The power spectrum of the reconstruction involves the four-point function of the observed CMB; it can be written in 
terms of the fully-reduced trispectrum as 



(CD = < +4^i j ^^g,,.,(L)5.3u(L) [T[-2jW +Tr[;:;:j(|li - I3I) +T^--j(|li - U\) 



(B17) 



where li I2 = L = I3 + I4. Here, ivf^ = Al arises from the Gaussi an (i.e. disconnected) part of the four-point 
function. To second order in (f>, the primary coupling term in Eq. (B17l gives the lensing power spectrum, Cf'^ that 
we aim to reconstruct. The other two couplings give the n'^^^ bias found in Ref. ED): 



(1) _ 



iV 



To fourth order in 
bias: 



Ili-UI 



ii • (ii - i4)(-i3) • (ii - i4)qt* , ,c««c?^ 



(B18) 



N 



(2) 



the primary coupling of the dominant trispectrum terms, Eqs. (B9| and (BIO I, give the iV|' 



(2) 



(2^) 



2 5l3l4(L)/l3l4(L) 



d'l 



5u,,(L) I j ^1 . L[l . (1 - h)rcrc;\^ Rh ■ LilcZ^ 



2 ^^^^^ 

This is the flat-sky version of Eq. |34|. We show in the text that ivf Vc^* 

can be expressed in terms of the leading- 
order difference between the lensed and unlensed power spectra. More correctly, the final term in Eq. (B19) can be 
shown to involve the difference between the power spectrum of the unlensed CMB and the cross spectrum of the 

lensed CMB with the lensed temperature gradient, but the latter is equal to the lensed spectrum to better than one 

(2) 

percent (Lewis, Challinor & Hanson, in prep.). We show further in the text that Nj^ ' /2 also arises when computing 

the power spectrum C^*^. Using techniques from Lewis, Challinor & Hanson (in prep.), the latter can be computed 
non-perturbatively in the lensing deflection as 



d% 
(27r)2 



51,1, (L) (li • LQ^^^ + I2 • LC7< 



(B20) 



where the approximation arises from replacing the exact cross power spectrum of the lensed temperature and the 

lensed temperature gradient with the exact Cf^. We expect {Cff') to equal Cff' plus half the non-perturbative 

generalization of 7V|^'' . We thus expect the dominant low-^ fractional bias in Cff to be controlled at higher order by 
the difference between the (exact) lensed spectrum and the unlensed spectrum. 
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